Baryon production and net-proton distributions in relativistic heavy ion collisions 
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The higher order moments of the net-baryon distributions in relativistic heavy ion colhsions are 
useful probes for the QCD critical point and fluctuations. We study the net-proton distributions and 
their moments in a simple model which considers the baryon stopping and pair production effects 
in the processes, ft is shown that a single emission source model can explain the experimental data 
well. Centrality and energy dependence of the distributions and higher moments is discussed. 
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I. INTRODUCTION 

The investigation of QCD phase diagram is fundamen- 
tal to our understanding of strong interactions. At van- 
ishing baryon chemical potential, lattice QCD calcula- 
tions predict the occurrence of a cross-over from hadronic 
phase to the quark-gluon plasma phase above a critical 
temperature of about 170-190 MeV [1, 2]. A distinct sin- 
gular feature of the phase diagram is the QCD critical 
point [3] which is located at the end of the transition 
boundary. A characteristic feature of the critical point is 
the divergence of the correlation length i^ and extremely 
large critical fluctuations. In ultra-relativistic heavy ion 
collisions, however, because of finite size and rapid ex- 
pansion of the system, those divergence may be washed 
out. As estimated in Q, the critical correlation length 
in heavy ion collisions is not divergent but about 2-3 fm. 
Remnants of those critical large fluctuations may become 
accessible in heavy ion collisions through an event-by- 
event analysis of fluctuations in various channels of con- 
servative hadron quantum numbers, for example, baryon 
number, electric charge, and strangeness [^. In an en- 
ergy scan there would be a non-monotonic behavior of 
non-Gaussian multiplicity fluctuations, which would be 
a clear signature for the existence of a critical point. In 
fact, at vanishing chemical potential it has been shown 
that moments of conservative charge fluctuations are sen- 
sitive indicators for the occurrence of a transition from 
hadronic to partonic matter J5|. 

Recently, the higher order moments of net-baryon dis- 
tributions in heavy ion collisions at RHIC energies have 
aroused great interest both experimentally [6[ and the- 
oretically [THlOj. Experimentally, neutrons can not be 
detected easily and the reconstruction efficiency is very 
low for strange hadrons. Fortunately, theoretical calcu- 
lations confirmed that the net-proton distribution can 
be a meaningful observable for the purpose of detect- 
ing the critical fluctuations of net baryons in heavy ion 
collisions [li|. The theoretical interest on these higher 
order moments comes from the discovery of the relation 
between the moments and the thermal fluctuations near 
the critical points. It is shown that the higher order 
moments have stronger dependence on the correlation 



length ^ and are therefore more sensitive to the critical 
fluctuations. If some memory of large correlation length 
persists in the thermal medium in hadronization process, 
this must be reflected in higher order moments of the dis- 
tributions. It has been predicted [13] that the third mo- 
ment, called skewness, is proportional to ^'*'^ and fourth 
moment, or kurtosis, proportional to ^^ while the sec- 
ond moment proportional to ^^. More importantly, the 
moments are closely related to the susceptibilities of the 
thermal medium. It has been argued that information 
of QCD phase diagram and the critical point can be ob- 
tained from the energy dependence of those moments [3] ■ 
The moments of net-proton distributions are studied with 
different theoretical models such as AMPT and UrQMD 
[9!], HIJING [lj|, and hadron resonance gas model [ij] 
etc. All those theoretical models are quite complicated 
and many microscopic processes are involved, and as a 
result many parameters can be tuned in the investiga- 
tions. Therefore, the underlying physics behind the ex- 
perimental results on the higher order moments of the 
net-proton distribution is not very transparent from the 
model studies. In addition, those studies focused on the 
moments only and made no direct comparison with the 
experimentally obtained distributions. 

In this paper, we will investigate the net-proton dis- 
tributions in Au-I-Au collisions at ^snn = 200GeV 
from very simple physics considerations: baryon stop- 
ping and baryon pair production. These physics effects 
are well known from studying heavy ion collisions in the 
past decades. We will show that such simple physics 
can be used to reproduce the experimentally observed 
net-proton distributions at different colliding centralities 
with parameters chosen properly. Then higher moments 
can be calculated numerically from the distributions. In 
this way the centrality dependence of those moments can 
be predicted. 

This paper is organized as follows. In next section, we 
will address the physics points in our considerations for 
an emission source. Analytical expressions for the net- 
proton distribution will be given. The model will be used 
to fit the experimental data on the distribution. The 
model can fit the data nearly perfectly. The centrality 
dependence of the moments will then be predicted. Also 



the moments at LHC energies are discussed. The last 
section will be for a brief summary. 



II. MODEL CONSIDERATION FOR AN 
EMISSION SOURCE 

It was well established that the net-baryon number 
would be zero in heavy ion collisions in central rapid- 
ity region if there were no nuclear stopping in the pro- 
cesses. Because the nuclear stopping effect depends on 
the collision energy and the size of the system, the baryon 
number stopped in a rapidity region is closely related to 
the number of participant nucleons A^part ■ In more cen- 
tral collisions the net baryon number will be larger. We 
consider a case in which all final state baryons are as- 
sumed being produced from one emission source. The 
initial mean nucleon number in the source is denoted as 
B which may be different for different colliding central- 
ities. Considering the randomness and independence of 
the nucleon-nucleon collisions in a heavy ion collision, the 
probability of finding A'o baryons stopped in the kine- 
matic region under investigation can be assumed, with 
the given mean number B, as 



Po(iVo,i3) = -— exp(-5) 

iVo! 



(1) 



Out of those A'o stopped nucleons, some of them are pro- 
ton, others neutrons. The probability of finding Np pro- 
tons from TVq nucleons is 



^P„W„ 
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(2) 



with p = Z/A the fraction of proton in the nucleus. The 
above formulas determine the distribution of net proton 
number in the source in the initial state of the collisions. 
Then one can consider the baryon production in the 
collisions. Baryons can be produced from various chan- 
nels. The baryon number is a conserved observable, thus 
baryons must be produced in baryon-antibaryon pairs. 
The pair production may be independent, and as a re- 
sult of the independence, the probability for producing M 
baryon pairs must be a Poissonian with the given mean 
number of produced pairs /i as a parameter 



M 



Po{M,^i) 






exp(-/i) 



(3) 



The parameter /i depends on the colliding central- 
ity. For more central collisions, the colliding system is 
larger, therefore fi should be larger, and more nucleon- 
antinucleon pairs can be produced in the process. 

It is the right place to compare the number distribu- 
tions used in this paper and in others. We use a Pois- 
son distribution for the baryon pair distribution, sup- 
posing the independent production of the pairs. In [16[ 
the distributions for both proton and anti-proton are 
assumed Poissonian, implying that protons and anti- 
protons are produced completely independently. There- 
fore the baryon number conservation may be violated in 



any event. In [T^l, a canonical ensemble is employed to 
derive the number distribution for tt systems. This is rea- 
sonable because there are a lot of tt particles in the final 
state of heavy ion collisions. But a simple transportation 
of the method to the case for baryon production may be 
problematic, because the relevant baryon particle num- 
ber may be not large enough for an equilibrium statistical 
description. 

In the strong production of nucleon-antinucleon pairs, 
isospin is conserved. Suppose that Ni protons, N2 anti- 
protons, N^ neutrons and 7V4 anti- neutrons are produced 
in the process, the conservation of isospin reads Ni — 
N2 — N3 — N4 if the effect from the presence of mesons 
is neglected. Thus we have Ni + N4 = N2 + N3 = M. 
The probability of finding A^i protons can be assumed as 



iiNi,M) 






(4) 



In writing this equality, we assume that all the produced 
pairs are within the kinematic range detected experimen- 
tally. Of course, this is a rather rough approximation. 
In fact, some of the produced nucleons can go out of 
that range and cannot be included in measuring the net- 
protons in the event. The effect from limited kinemati- 
cal acceptance can be taken into account by introducing 
one more parameter for the probability of the produced 
baryon in the detected region. To avoid this complexity, 
in this paper, such an effect is effectively treated as hav- 
ing the number of pairs M a little smaller in the event. 
Therefore, the value of the parameter /x obtained from 
the fitting in this paper should be a little bit smaller 
than the real one. In the same way, the probability of 
finding iV2 anti-proton is Qi{N2,M). Then the distri- 
bution of net-proton Ap from an emission source can be 
expressed as 



P{Ap) 



No,Np,M,Ni,N2 



PoiNo,B)QoiNo,Np) 



Qi{Ni,M)Qi{N2,M)Po{M,fi)SAp,N,+N^-N, ■ (5) 

By inserting an identity expression Sm^n = 
Jp'^ dxe*^™~"^^/27r, the above equation can be rewritten 



P{Ap) = / — 
Jo I" 



-{2Bp+)j.) sin^ 



2 cos(xAp — Bp sin x){6) 



As can be seen from the above expression, the net-baryon 
distribution depends on two combined parameters, Bp 
and p instead oi B, p and p separately. One can check 
easily that Bp is the mean value of the distribution 
PiAp). 



III. COMPARISON WITH THE 
EXPERIMENTAL DATA 

The expression Eq. ^ enables us to compare the cal- 
culated net-proton distributions from an emission source 
to the experimental data from STAR 6|, as shown in Fig. 
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Bp 


0-5% 


351.4 


14.5 


1.65 


30-40% 


114.2 


5.5 


0.632 


70-80% 


13.4 


0.83 


0.075 



TABLE I. Fitted parameters for Fig. [T] 



[TJ The parameters used are tabulated in TABLE ID The 
parameters show the expected behaviors from central to 
peripheral collisions. As one can see from the figure, the 
agreement with the data is very good over five orders of 
magnitude. Almost all calculated points for F{Ap) lie 
within the experimental error bars. 
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FIG. 1. The net-proton distributions from single emission 
source. The points are from Ref. 61, and the curves are 
calculated from Eq. |6] 



FIG. 2. The mean net-proton from multiple emission sources. 
The points are from Ref. [y|, and the curve is from our model 
calculation. 
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To make predictions for the net-proton distributions at 
other centralities, one can parameterize the values of pa- 
rameters ^ and Bp tabulated in TABLE|T]by polynomials 
of the number of participants Np as 



H = 0.171 + 0.0495iVp - 2.5 x 10"^iV„^ 



(7) 
Bp = (-4.63 + 5.99A^p - 3.7 x 10"3)/1000. (8) 



From this parameterization, one can calculate the mo- 
ments for the distributions easily. The obtained moments 
are shown, as functions of the number of participants 
-^part, in Figs. [2][5l The corresponding experimental data 
from ^] are shown in the figures for comparison. The cal- 
culated mean, variance, skewness and kurtosis are well 
in agreement with the data. The good agreement shows 
that the basic merits for the baryon production mecha- 
nism have been exhibited in our model consideration. 

In Ref. [ifl], the observed moments are related to the 
ones from one emission source by using the central limit 
theorem. This expectation is deduced from the indepen- 
dent emission of particles from each source. In fact, if 
the baryons are produced from Ns identical sources, one 



FIG. 3. The variance for the net-proton distributions from 
multiple emission sources. The points are from Ref. [g, and 
the curve is from our model calculation. 
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FIG. 4. The skewness for the net-proton distributions from 
multiple emission sources. The points are from Ref. [g, and 
the curve is from our model calculation. 
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FIG. 5. The kurtosis for the net-proton distributions from 
multiple emission sources. The points are from Ref. [g|, and 
the curve is from our model calculation. 



can get relations of the moments for the measured distri- 
butions and those for the emission sources as 

M = A'hNs, 
a ^ Oi \fN~s, 
S - SJy/N^, 

K = Ki/Ns, 

where quantities with subscript i are for moments from 
one emission source. From the above expressions on cen- 
trality dependence, one can expect constant Sa and na"^ 
for all iVpart- In our fitting with a single emission source, 
both n changes strongly with centrality. As a result of 
the centrality dependence of Bp and /i, the centrality 
dependence of the moments are well reproduced, as can 
be seen from Figs. [2][5l One can also calculate moment 
products Sa and kct^ . The centrality dependence of the 
products are shown in Fig. [5] In the A^part range shown, 
Sa increases a few percent, while Ka^ is almost ex actly 
1, as expected from the hadron resonance gas model jlSj . 

With the good agreement with STAR data at RHIC 
energy at hand, one can go one step further to pre- 
dict the higher order moments of the net-proton distri- 
butions at LHC energies. For specific, let the center- 
of-mass energy of the colliding nucleon-nucleon pair be 
2.76TeV. From fitting data in [6[ one can find the center- 
of-mass energy dependence of the parameter Bp and ex- 
trapolate to LHC/ALICE energy. In this way, one gets 
Bp ~ at y/sjMN = 2.76 TeV. To get the parameter p at 
LHC/ALICE, one can write p ex exp(?Tip/r), with T be- 
ing the effective emission temperature of the sources and 
nip the mass of proton. In [18| the effective temperature 
of the medium is given as a function of the center-of-mass 
energy y^ as 
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FIG. 6. The moment products, Sa and Ka^, for the net- 
proton distributions from multiple emission sources as func- 
tions of A'part. The dotted line is for kct^ = 1 expected from 
the hadron resonance gas model. 



with TLim — 164 MeV. At LHC, the center of mass en- 
ergy is much higher than at RHIC, so the value of pa- 
rameter p is much smaller than obtained in the above. 
The smaller value of p will allow more baryon pairs to 
be produced from a source. The predicted moments are 
shown in Fig. [71 Now the mean value of the distribution 
is zero. Because of zero mean value, the net-proton dis- 
tribution at LHC/ALICE is symmetric about and one 
gets zero value for the skewness. Since the kurtosis is 
zero for Gaussian distributions, extremely small k value 
at large A'^part at LHC energies means that the net-proton 
distributions at LHC energies can be well parameterized 
by Gaussian, but not for small A^part- 
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FIG. 7. The moments for the net-proton distribution from 
multi-emission sources at ^snn = 2.76 TeV. 
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,/[l + exp(2.6-ln(v^)/0.45)] 



IV. CONCLUSION 



The higher order moments of the net-proton distri- 
butions in relativistic Au-|-Au colhsions at ^/sNN = 
200GeV are studied from a simple model with effects 
from initial baryon stopping and final baryon pair emis- 
sion taken into account. We have demonstrated that by 
employing a single emission source model, the distribu- 
tions at different collision centralities can be well repro- 
duced. Then the higher order moments for the distribu- 
tions can be calculated without new free parameters. The 
calculated moments agree well with the experimental re- 
sults. The predicted moments for LHC Pb+Pb collisions 
need to be verified experimentally. 

It should be mentioned that nothing else is assumed in 
this model except an initial net-proton and a finite prob- 



ability for producing baryon pairs from sources. There- 
fore, our model has nothing to do with thermal equi- 
librium and/or critical fluctuations. Because our model 
consideration is based on normal physics effects, our re- 
sults can be used as a baseline for detecting novel physics 
in the processes. 
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